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ABSTRACT 

In  this  paper  we  consider  the  rational  canonical  form  of  arbitrary 
polyphase  matrices  and  use  it  to  derive  a  simple  implementation  of 
paraunitary  filter  banks  (PUFBs)  based  on  a  cascade  of  elementary 
building  blocks.  Furthermore,  this  decomposition  is  shown  to  be 
easily  extendable  to  include  a  large  class  of  perfect  reconstruction 
filter  banks  (PRFBs)  and  can  be  especially  useful  for  deriving  the 
initial  condition  of  PUFB  design  algorithms. 

Index  Terms —  rational  canonical  form,  polyphase  matrices, 
paraunitary  filter  banks 

1.  INTRODUCTION 

The  polyphase  matrix  of  an  M-channel  causal  finite  impulse 
response  (FIR)  filter  bank  may  be  written  as 

H(z)  = 

k= 0 

where  z  =  re JC0  and  Hk  is  a  scalar  M-by-M  matrix.  The  matrix  H(z) 
resides  within  the  ring  of  M-by-M  matrices  with  Laurent 
polynomial  entries,  denoted  by  M (N,  C [z,  z *]),  with  C  the  field  of 
complex  numbers.  The  rational  canonical  form  of  this  matrix  is 
not  guaranteed  to  exist  within  the  Laurent  polynomial  ring  C[z,z4] 
since  this  ring  is  not  a  field.  In  this  paper  we  will  show  that  the 
rational  canonical  form  of  H(z)  does  indeed  exist  as  another  matrix 
in  the  ring  M(JV,  C[z,  z  1).  If  H(z)  corresponds  to  the  polyphase 
matrix  of  a  paraunitary  filter  bank  (PUFB),  then  the  algorithm  for 
computing  the  rational  canonical  form  illustrated  in  this  paper  may 
be  used  to  implement  H(z)  as  a  particularly  simple  product  of 
elementary  building  blocks.  Furthermore,  this  new  decomposition 
of  H(z)  may  be  used  to  represent  a  large  class  of  perfect 
reconstruction  filter  banks  (PRFBs)  as  well  and  can  be  used  to 
provide  initial  conditions  for  PUFB  design  algorithms. 

2.  EXISTENCE  OF  RATIONAL  CANONICAL  FORM 

The  rational  canonical  form  of  a  polyphase  matrix  H(z)  over  the 
ring  C[z,  z'1]  exists  as  a  matrix  over  the  field  K.  K  denotes  the 
field  of  rational  functions  p(z)/q(z),  q(z)f 0,  with  coefficients  in  C. 
The  field  K  is  the  field  of  fractions  of  the  ring  C[z,  z"1].  This  field 
is  also  denoted  C(z).  Any  matrix  H(z)  in  M(iV,  C[z,  z"  ])  is  similar 
to  a  matrix  in  M (N,  K)  in  rational  canonical  form.  In  other  words, 
there  exists  an  invertible  M-by-M  matrix  P(z)  such  that  P(z) 
'H(z)P(z)  is  in  rational  canonical  form.  A  matrix  R(z)  in  rational 
canonical  form  is  a  matrix  such  as 

Cj(z)  0  ■■■  0 


0  0  0  CM  (z) 


where 


C,(z)  = 


0  0 
1  0 


0  0 


~Po(z)/4o(z) 

~  Pi(z)/ti(z) 


e  M(A,K). 


The  matrix  Cj(z)  is  the  companion  matrix  corresponding  to  the 
invariant  factor  afx)  of  H(z).  The  invariant  factors  a;(x)  may  be 
made  unique  by  requiring  that  pfz)  and  qfz)  are  monic 
polynomials. 

The  characteristic  polynomial  of  H(z)  is  a  polynomial  in  the 
indeterminant  x  with  coefficients  in  the  Laurent  polynomial  ring 
C[z,  z"1],  namely  c(x)  =  det[xl  -  H(z)].  The  invariant  factors  of 
R(z)  are  the  monic  factors  of  c{x)  [  1 J .  The  invariant  factors  of 
H(z)  and  R(z)  are  the  same,  namely  the  monic  polynomials  afx). 

The  important  contribution  of  this  section  is  the  following 
claim. 

Claim :  The  rational  canonical  form  of  the  polyphase  matrix  H(z) 
is  a  matrix  with  entries  in  the  ring  C[z,  z  '].  In  other  words,  it  is  an 
element  in  the  ring  M(A,  C[z,  z"1]). 

Proof.  Since  C  is  a  field,  the  ring  C[z,  z1]  is  a  unique  factorization 
domain  (UFD).  The  field  K  =  C(z)  is  the  field  of  fractions  of 
C[z,  z  1 1 .  The  characteristic  polynomial  c(x)  is  in  the  ring 
C[z,  zqjjc].  By  Gauss’  Lemma,  since  c(x)  can  be  factored  in 
C(z)M,  then  it  is  reducible  in  C[z,  z"1] [x].  Consequently,  the 
invariant  factors  of  a  matrix  H(z)  in  M (N,  C[z,  z"1])  are 
polynomials  in  the  ring  C[z,  z"1]  and  therefore  the  rational 
canonical  form  R(z)  is  a  matrix  in  the  ring  M(N,  C[z,  z4]). 


3.  AN  ALTERNATIVE  FORM  FOR  R(z) 

Claim :  The  rational  canonical  form  R(z)  of  a  polyphase  matrix 
H(z)  can  be  written  as  a  Laurent  polynomial  with  matrix 
coefficients.  In  other  words, 

"C,(z)  0  ■■■  0 

R(z).  °  ^  /  1 

_  0  0  0  CM(z)_ 

=  R„  +  R,z4  +  R2z-2  +  ■  ■  ■  +  Rnz~w. 

Proof.  The  claim  follows  once  one  proves  that  the  ring  of  matrices 
with  Laurent  polynomial  entries  is  isomorphic  to  the  ring  of 
Laurent  polynomials  with  matrix  coefficients.  It  is  sufficient  to 
prove  that  [2], 

M(v,C[jc])  =  M(v,  C)[je] 

which  implies  that 

m(n,  c[z,z4D  =  M(N,  C)[z,z4]  . 

The  proof  is  omitted  for  brevity. 
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4.  COMPUTATION 

To  compute  the  rational  canonical  form  R(z)  of  the  M-by-M  matrix 
H(z),  one  approach  is  to  diagonalize  the  matrix  xl  -  H(z)  [1],  The 
invariant  factors  of  H(z)  will  then  appear  on  the  diagonal.  The 
invariant  factors  should  be  monic  and  satisfy  the  divisibility 
condition  a;(x)  divides  ai+1(x).  The  direct  sum  of  the  companion 
matrices  associated  with  each  invariant  factor  yields  the  matrix 
R(z)  in  rational  canonical  form.  By  keeping  track  of  the  row 
operations  used  to  diagonalize  xl  -  H(z)  one  can  also  construct  the 
matrix  P(z)  such  that  P(zy‘H(z)P(z)  =  R(z). 

The  following  three  elementary  row  and  column  operations 
can  be  used  to  diagonalize  xl  -  H(z): 

1.  interchange  two  rows  or  columns 

2.  add  a  multiple  in  K[x]  of  one  row  or  column  to  another, 
e.g.  add  p(x)  times  the  /'th  row  to  the  /th  row 

3.  multiply  any  row  or  column  by  a  unit  in  K[x],  i.e.  by  a 
nonzero  element  in  K. 

The  matrix  P(z)  can  also  be  computed  systematically  [1],  First,  let 
dy...,dm  denote  the  degrees  of  the  monic  nonconstant  polynomials 
ai(x),...,am(x)  appearing  on  the  diagonal.  Begin  with  the  matrix 
S  =  I,  the  identity  matrix.  For  each  row  operation  used  to 
diagonalize  xl  -  H(z),  change  the  matrix  S  as  follows. 

1.  If  the  /th  and  /th  rows  were  interchanged,  then 
interchange  the  /th  and  jth  columns  of  S, 

2.  if  Row,  +  p(x)Row]  A  Row-,,  then  subtract  the  product  of 
the  matrix  p(H(z))  times  the  /th  column  of  S  from  the 
y'th  column  of  S,  i.e.  Co/j  -  p(H(z))Co/;  Co/j, 

3.  if  the  /th  row  is  multiplied  by  a  unit,  u,  then  divide  the 
/th  column  of  S  by  u. 

Once  the  matrix  xl  -  H(z)  has  been  diagonalized,  the  first  M  -  m 
columns  of  S  will  be  zero.  Then  for  each  i  =  1 multiply  the 
/th  nonzero  column  of  S  successively  by  H(z)°  =  I.  H(z),  H(z)2,..., 
Hfz)*’1,  where  d-,  is  the  degree  of  a;(x).  Use  the  resulting  column 
vectors  in  this  order  as  the  next  d-,  columns  of  a  matrix  P(z).  Then 
P(z)-‘H(z)P(z)  =  R(z). 

5.  APPLICATION  TO  PUFBS 

In  this  section  we  consider  the  example  of  paraunitary  filter  banks. 
The  analysis  (or  synthesis)  polyphase  matrix  of  a  normalized 
PUFB  satisfies  the  condition, 

H(z)h(z)"  =  I  for  z  =  e1". 

Furthermore,  M-channel  PUFBs  may  be  decomposed  into  a 
product  of  elementary  building  blocks  as  in  [3], 

H(®)  =  GtA(®)---G,A(o>)QJ  (1) 

where  L  is  the  Smith-McMillan  degree  of  H(co),  Gk  and  Q  are 
M-by-M  orthogonal  matrices,  A(co)  =  diag(I,  for  M  even,  and 
J  =  diag(±l,...,±l).  Each  matrix  Gk  and  Q  can  be  written  as  the 
product  of  VzM(M  -  1)  Givens  rotation  matrices  in  the  sequence, 

=  {b„_2,„-i}"Km-1"B12Km-1-B01}.  (2) 

The  matrix  B,,  corresponds  to  a  Givens  rotation  matrix  with  cos(0n) 
in  the  /th  row  and  /th  column,  sin(0n)  in  position  (/,  j),  -sin(0„)  in 
position  (j,  /),  and  cos((9n)  in  position  (j,j)  with  1  <  n  <  ZiM(M  -  1). 


5.1.  CaseZ  =  1 ,  M  -  2 

For  a  simple  example,  let  M  =  2  and  L  =  1  in  (1)  above.  Without 
loss  of  generality  set  QJ  =  I.  Then  the  polyphase  matrix  is 


H(z)  =  G,A(z)  = 


cos(#) 

-sin(0) 


sin(<9)ir  1 
cos(0)_|_O 


(3) 


Table  1  lists  the  steps  to  compute  the  rational  canonical  form  of 
H(z)  and  the  change  of  basis  matrix  P(z).  In  this  table,  R,  denotes 
the  /th  row,  and  Cj  denotes  the  /th  column. 


Table  1.  Sample  Rational  Canonical  Form  Computations 


Steps  to  compute  R(z) 

Steps  to  compute  P(z) 

Step  1:  Form  the  matrix 
xl  -  H(z) 

Step  1:  Form  the  matrix 

S  =  I 

Step  2:  — sin(f?)i?!  A  Ru 
(x  -  cos (6))R2  R2 

Step  2:  -(l/sin(0))Ci ->  Ci, 

(H(z)  -  cosC^UC,  +  C2AC, 

Step  3:  -Ri  +  R2~)  R2 

Step  3:  sin(0)Ci  ->  Ci 

Step  4:  (l/sin(0))f?[  R\ 

Step  4:  (z/sin {0))C2  ->  C2, 

-Ci  ^  Ci 

Step  5:  z  lsin(8)R2  ->  R2, 

(x2 - xcos(0)(l+z"‘)  +  z')Ri 

Ri 

Step  5:  Ci  -  (H(z)2  - 
cos(60H(z)(l  +  z ')  +  z"‘l)C2 

c2 

Step  6:  -Rt  Ri 

Step  6:  (H(z)  -  cos(0)I)Ci  -> 

Cu  sin  (6»)z  ‘C2  A  C2 

Step  7:  Ri  +  R2  Ri 

Step  7:  (H(z)2-cos(6»)H(z)(l + 
z')  +  zll)Ci  C, 

Step  8:  (l/(x  -  cos(0))f?i  -A  Rlt 
(l/(z'sin(6»))«2  R2 

Step  8:  Now, 

[0  0  1 
s  = 

0  -2z~‘  sin(6>)2 

Also,  m  =  1,  dt  =  2. 

Step  9:  (l/(x2  -  xcos(<9)(l  +  z"1) 

+  z')Ri^Rt 

Step  9:  Form  the  matrix, 
P(z)=[lC2  H(z)C2] 

The  result  of  these  computations 

„n-r°  _z_1 

|_1  cos(0)(l  +  z-‘) 

S 

0  0  1  Jo  -1 

z  1 

1  cos(^)  '  0  cos(&) 

P(z)  =  -2sin(0):z-1 


z  ‘sin(0) 
Z-1  cosffl) 


(4) 


V(zY  =-y2zcsc(e)2 


-cot(^) 

zcsc(d) 


1 

0 


Furthermore,  P(zy‘H(z)P(z)  =  R(z). 

Using  the  rational  canonical  form  of  H(z),  one  can  write 
H(z)  =  H,(z)  +  z~‘H2(z)  where  Hi(z)  and  H2(z)  are  rank  one 
matrices  defined  as 


H,  z  = 


H2(z)  = 


«(*) 


cos! 
cos(#)co  t(0) 
cos(#) 

CSC  (#)(cos(#y 


z  tan! 


(4 


[1  0], 


Now  for  M  =  2,  the  decomposition  of  PUFBs  given  in  (1)  can  be 
rewritten  as, 


HW  =  n(HliW  +  z-‘H2i(z)). 
1=1 


6.  SIMPLE  IMPLEMENTATION  OF  PUFBS 

Using  the  decomposition  Pi(z)Ri(z)Pj(z) 1  =  G;A(z)  for  each 
building  block  of  a  2-by-2  polyphase  matrix  in  (1),  a  two  channel 
PUFB  may  be  constructed  as  a  cascade  of  the  lattice  sections 
shown  in  Fig.  1.  Here  Vi,  v2  are  the  input  signals  and  u\,  u2  are  the 
output  signals.  This  lattice  structure  lends  itself  easily  to  an 
implementation  in  hardware. 
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Fig.  1.  Lattice  Implementation  of  PUFB 
Similar  lattice  structures  can  be  derived  for  PUFBs  with 
greater  than  M  =  2  channels  using  (2),  the  rational  canonical  form 
for  (3)  and  the  rational  canonical  form  for  a  scalar  Givens  rotation 
matrix.  Using  the  same  computational  algorithm  described  in 
Section  4,  a  scalar  Givens  rotation  matrix  can  be  written  in  rational 
canonical  form  with  PRP  1  =  G  as, 

0  —1  cos 

1  2cos($)  '  -sin(P)  cos(P) 


R  = 


P  1  = 


2cos((?)  ’ 
-  cot(<?)  1 

esc  fa)  0 


.  P  = 


0  sinl 
1 


in(<?)' 

(*). 


(5) 


6.1.  Case  M  -  3 

Using  (2),  a  3-by-3  paraunitary  polyphase  matrix  may  be 
written  as. 


h(z)  =  B12  fa  )bo,  fa2  )boi  fa3  )a(z), 

with  A (z)  =  diag(l,  z  l,  z "').  A  lattice  implementation  can  be 
constructed  from  the  following  matrices  derived  using  (4)  and  (5), 

0  o' 


cscfa )  0 

-cotfa2)  0  1 

0  1  0 
cscfa2 )  0  0 


”1  0 

0 

“l  0 

o  in 

Bufa)  = 

0 

D 

sinfe) 

0 

0 

-1  0 

0 

cos(^j ) 

0 

1 

2cos(^j)  0 

0 

0 

sin(^2 ) 

r° 

0 

-1 

B„2(ft)  = 

0 

1 

0 

0 

1 

0 

1 

0 

COs(^2  ) 

1 

0 

2  cos  (Or,) 

0 

sin(^3) 

o' 

'1 

0  0  1 

®0lfe)^(z)_ 

1 

cos(6?3, 

0 

0 

z_1  0 

0 

0 

1 

0 

0  z_1  J 

0  -  z~l 

1  cnslr^-  )(l  +  z 
0  0 


—  cot(<% )  1  O' 
zcscfa3)  0  0 
0  0  1 


Similar  matrices  and  lattice  structures  may  be  derived  using  (2), 
(4),  and  (5)  for  paraunitary  polyphase  matrices  with  M  =  4,  or 
larger.  Those  examples  will  not  be  presented  here  to  conserve 
space. 


7.  PERFECT  RECONSTRUCTION  FILTER  BANKS 

Consider  the  case  M  =  2  again.  A  perfect  reconstruction  filter 
bank  satisfies  the  property  that  H(z)E(z)  =  czm I,  where  H(z)  is  the 
analysis  polyphase  matrix  and  E(z)  is  the  synthesis  polyphase 
matrix,  c  f  0  and  m  is  an  integer.  A  PUFB  is  the  special  case 
where  E(z)  =  Hh(z).  A  perfect  reconstruction  filter  bank  may  be 
constructed  by  parameterizing  the  rational  canonical  form  R(z)  and 
the  conjugation  matrix  P(z)  independently  using  one  angle  for 
each.  For  example,  the  system  with 


H(z) 


0  sin(0, ) 

"0  -1 

-cos(Pj)  sin(0,) 

‘l  0  ' 

1  cos(P, ) 

1  2cos(#2) 

1  0 

_°  z_1. 

CSC 


fa) 


e(z)  = 


'1  o' 

0  sin(0,) 

2cos(#2)  1 

0  z 

1  cos(P, ) 

-1  oj 

csc(#j)  (6) 


-cos  fa)  sin  faJl 
1  0 

is  PR  since  E(z)  =  H(z)’1.  Note  that  for  the  special  case  when  8\  = 
62 ,  H(z)  is  paraunitary.  Therefore,  this  set  of  PRFBs  constructed 


using  two  parameters  for  every  building  block  includes  the  set  of 
PUFBs  as  a  proper  subset. 

8.  APPLICATION  TO  PUFB  DESIGN 

Since  the  space  of  PRFBs  is  larger  than  and  includes  the  space  of 
PUFBs,  one  can  design  a  PUFB  by  searching  over  the  space  of 
PRFBs  for  an  optimal  solution  and  then  choosing  the  paraunitary 
filter  bank  closest  to  it  for  an  initial  condition  in  a  strictly  PUFB 
search.  The  advantage  of  this  approach  is  that  if  the  initial 
condition  of  a  strictly  PUFB  search  is  closer  to  the  optimal 
solution,  then  the  more  likely  it  is  that  a  nonlinear  optimization 
program  will  settle  on  the  global  solution  instead  of  a  local 
minimum.  Consider  the  following  PUFB  design  algorithm  which 
minimizes  the  Mean  Square  Error  (MSE),  £  between  a  desired 
polyphase  matrix,  D(ro),  and  an  approximation  HRu). 

8.1.  PUFB  Design  Algorithm 

Step  1.  Form  the  MSE  objective  function, 

2  K 

mse(dM,p(®))  =  <r  J  wVtoto-Xvld* 

0 

where  II .  IIF  denotes  the  Frobenius  norm,  W(co)  is  a  scalar 
weighting  function  set  equal  to  one,  P(co)  is  an  M-by-M  perfect 
reconstruction  polyphase  matrix  as  described  in  Section  7,  and 
D(ro)  is  the  desired  or  ideal  polyphase  matrix. 

Step  2.  With  the  angles  0\  as  free  parameters,  minimize  < f  over  all 
PRo)  using  a  nonlinear  optimization  program  such  as  the  Broyden- 
Fletcher-Goldfarb-Shanno  (BFGS)  algorithm.  Call  the  optimal 
solution,  P*Ro). 

Step  3.  Compute  the  paraunitary  polyphase  matrix  H’Ro), 
parameterized  as  in  (1)  and  (2),  which  minimizes 
MSE(P*Ro),  HRo))  over  all  HRo). 

Step  4.  Use  H'(co)  as  the  initial  starting  point  in  a  nonlinear 
optimization  program  which  minimizes  MSE(D(co),  HRo))  over  all 
HRo).  Call  this  optimal  solution  H*Ro).  Denote 
MSE(DRo),  H*Ro))  by  MSE0. 

The  advantage  of  using  the  proposed  PUFB  design  algorithm  is 
that  it  satisfies  the  following  property. 

Claim:  Define  MSE[  =  MSE(D(ro),  H'(co)).  Using  the  above  filter 
bank  design  algorithm, 

MSE,  <  MSE0  +  e, 

MSe(h'(<w).  H  ■  (<»))  <  2  ■  MSE0  +  £ 

where  c  is  a  small  nonnegative  constant. 

Proof:  Since  the  MSE  is  a  norm  on  the  space  of  polyphase 
matrices,  by  the  triangle  inequality, 

2/r  2k 

£  }  ||D(ffl)  -  H'H|  \dco  <^J  ||d(®)  -  P»|£  dco+ 

0  0 

2k 

0 

Since  P*(co)  is  the  optimal  solution  over  a  larger  set  of  matrices 
than  H*(co), 

2jt  In 

|  |d(o)  —  P’  da)  <  A-  j"  ||d(o)  —  H*  (ta)||f  da). 

0  0 

As  a  result, 
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2  n  In 

it  <  -h  da+ 

o  o 

2n 

1tJ||p!M-H'(®)||/®, 

o 

or  in  other  words,  MSE]  <  MSE0  +  s.  The  proof  that 
MSE(H'(co),  H*(to))  <  2-MSE0  +  e  is  similar. 

8.2.  Case  M  =  4 

For  M  =  4,  a  PR  polyphase  matrix,  P(z),  which  can  be  used  to 
design  a  PUFB  using  the  previously  described  algorithm  is 
constructed  as, 

P(z)  =  G1(ei,...,012)A(z)  = 

T23  (< 9 , ,  02  )t,  3  (03 ,  04  )t12  (05 , 06  )Tq3  (07 ,  08  )T02  (0, ,  0IO  )T01  (0, , ,  012  )a(z  ). 

The  inverse  of  P(z)  is  given  by, 

P(zt1  =  A(zfBol(0ll,012)Bo2(09,0IO)Bo,(07,08)BI2(0s,06)B1,(03,04)B23(01,02) 
with  the  matrices  By  as  in  (5)  and  (2)  except  with  two  angle 
parameters  as  in  (6).  The  system  so  described  is  PR  since 
P(z)P(z)  1  =  I.  The  matrix  T23  is, 


"l 

0 

0 

0 

T 

0 

0 

O' 

"l 

0 

0 

o' 

0 

i 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

0 

0 

sinfa) 

0 

0 

2cos(07) 

1 

0 

0 

-  cot(0, ) 

1 

0 

0 

l 

cos(0, ) 

0 

0 

-1 

0 

0 

0 

csc(0,) 

0 

The  matrices  T13,  T12,  T03,  T02,  Tra  are  constructed  similarly  using 
two  angle  parameters  for  each.  Note  that  if  =  for  i  odd,  then 
P(z)  is  paraunitary. 

9.  RESULTS 

The  PUFB  design  algorithm  was  used  to  design  a  finite  impulse 
response  PU  approximation  to  an  ideal  principal  component  filter 
bank  (PCFB).  PCFBs  are  described  in  detail  in  [4],  [5],  and  [6],  A 
PCFB  is  the  solution  to  the  problem  of  finding  optimal  Q-by-M 
and  M-by-Q  analysis  and  synthesis  polyphase  matrices,  with 
Q  <  M,  such  that  the  time-averaged  mean  squared  error  between 
the  vector  input  to  the  filter  bank  and  the  vector  output  is 
minimized.  An  ideal  PCFB  corresponding  to  an  infinite  order  PU 
filter  bank  has  channel  filters  with  brick-wall  responses.  This  ideal 
filter  bank  can  be  approximated  using  a  FIR  PU  filter  bank  by 
minimizing  the  mean  squared  Frobenius  norm  error  between  the 
desired  polyphase  matrix  of  the  ideal  PCFB,  D(o>),  and  the  FIR  PU 
synthesis  polyphase  matrix,  H(co).  PCFBs  are  an  interesting 
design  example  because  they  are  optimal  for  maximizing  coding 
gain  and  minimizing  mean-squared  error  in  the  presence  of 
quantization  noise.  It  has  been  proven  that  they  are  also  optimal 
for  any  concave  function  of  the  subband  variance  vector  [6]. 

Table  2  lists  the  performance  of  the  PUFB  design  algorithm 
and  compares  it  to  another  elegant  design  algorithm  described  by 
Tkacenko  in  [4],  Tkacenko  generously  made  available  the 
MATLAB  code  to  duplicate  his  results  on  the  Internet.  The  PR 
polyphase  matrix  was  populated  with  incrementally  more  two- 
parameter  building  blocks.  As  the  table  shows,  the  MSE  of  the 
optimal  PR  solution  decreased  monotonically  with  greater  degrees 
of  freedom,  until  the  number  of  free  parameters  was  too  great  for 
the  algorithm  to  converge.  A  PU  approximation  to  the  optimal  PR 
solution  with  the  largest  number  of  free  parameters  was  then  used 
as  the  initial  condition  for  a  search  over  PU  matrices  to  arrive  at 
the  final  filter  bank  solution.  The  length  of  the  channel  filters  in 
the  final  PUFB  was  eight  taps. 


Table  2.  Performance  of  PUFB  Design  Algorithm 


Algorithm 

MSE(D(ra),  H*(co)) 

Tkacenko  PU  solution 

1.7024 

PU  search  with  random  start  value 

1.5747 

PR  matrix  with  14  angles 

1.5578 

PR  matrix  with  16  angles 

1.5549 

PR  matrix  with  1 8  angles 

1.5539 

PR  matrix  with  20  angles 

1.5415 

PR  matrix  with  22  angles 

1.5395 

PR  matrix  with  24  angles 

3.6687  (algorithm  did 
not  converge) 

PU  search  with  PU  initial  condition 

1.5741 

As  Table  2  shows,  the  lowest  MSE  for  a  PUFB  was  attained  with 
the  PU  initial  condition.  Figure  2  illustrates  the  frequency 
response  of  one  of  the  optimal  channel  filters  (red),  and  overlays 
the  ideal  filter  response  (blue  dotted  line),  with  the  Tkacenko 
solution  (black). 


Fig  2.  Channel  Filter  One 

As  the  curves  show,  the  derived  filter  is  an  excellent 
approximation  to  the  ideal  brick-wall  filter  with  high  frequency 
selectivity  in  the  narrow  bandpass  regions. 

10.  CONCLUSIONS 

In  this  paper  we  showed  that  the  rational  canonical  form  of  a 
polyphase  matrix  exists  in  the  ring  M (N,  C[z,  z"1])  and  can  be 
written  as  a  matrix  in  the  ring  M(/V,  C)[z,  z"1].  The  factorization 
of  paraunitary  polyphase  matrices  into  a  similarity  transformation 
of  the  rational  canonical  form  makes  possible  a  very  simple  lattice 
implementation  of  the  filter  bank.  Also,  this  decomposition  lends 
itself  easily  to  a  representation  of  a  large  class  of  perfect 
reconstruction  filter  banks  and  can  be  used  to  find  an  initial 
condition  for  PUFB  design  algorithms  that  ultimately  yields  better 
MSE  performance. 
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